Predictive Power of Molecular Dynamics Receptor Structures in Virtual Screening

Molecular dynamics (MD) simulation is a well-established method for understanding protein dynamics. Conformations from unrestrained MD simulations have yet to be assessed for blind virtual screening (VS) by docking. This study presents a critical analysis of the predictive power of MD snapshots to this regard, evaluating two well-characterized systems of varying flexibility in ligand-bound and unbound configurations. Results from such VS predictions are discussed with respect to experimentally determined structures. In all cases, MD simulations provide snapshots that improve VS predictive power over known crystal structures, possibly due to sampling more relevant receptor conformations. Additionally, MD can move conformations previously not amenable to docking into the predictive range

In silico assessment of potential druggable pockets on the surface of α1-Antitrypsin conformers

The search for druggable pockets on the surface of a protein is often performed on a single conformer, treated as a rigid body. Transient druggable pockets may be missed in this approach. Here, we describe a methodology for systematic in silico analysis of surface clefts across multiple conformers of the metastable protein α1-antitrypsin (A1AT). Pathological mutations disturb the conformational landscape of A1AT, triggering polymerisation that leads to emphysema and hepatic cirrhosis. Computational screens for small molecule inhibitors of polymerisation have generally focused on one major druggable site visible in all crystal structures of native A1AT. In an alternative approach, we scan all surface clefts observed in crystal structures of A1AT and in 100 computationally produced conformers, mimicking the native solution ensemble. We assess the persistence, variability and druggability of these pockets. Finally, we employ molecular docking using publicly available libraries of small molecules to explore scaffold preferences for each site. Our approach identifies a number of novel target sites for drug design. In particular one transient site shows favourable characteristics for druggability due to high enclosure and hydrophobicity. Hits against this and other druggable sites achieve docking scores corresponding to a Kd in the µM–nM range, comparing favourably with a recently identified promising lead. Preliminary ThermoFluor studies support the docking predictions. In conclusion, our strategy shows considerable promise compared with the conventional single pocket/single conformer approach to in silico screening. Our best-scoring ligands warrant further experimental investigation

A combinatorial TIR1/AFB–Aux/IAA co-receptor system for differential sensing of auxin

The plant hormone auxin regulates virtually every aspect of plant growth and development. Auxin acts by binding the F-box protein transport inhibitor response 1 (TIR1) and promotes the degradation of the AUXIN/INDOLE-3-ACETIC ACID (Aux/IAA) transcriptional repressors. Here we show that efficient auxin binding requires assembly of an auxin co-receptor complex consisting of TIR1 and an Aux/IAA protein. Heterologous experiments in yeast and quantitative IAA binding assays using purified proteins showed that different combinations of TIR1 and Aux/IAA proteins form co-receptor complexes with a wide range of auxin-binding affinities. Auxin affinity seems to be largely determined by the Aux/IAA. As there are 6 TIR1/AUXIN SIGNALING F-BOX proteins (AFBs) and 29 Aux/IAA proteins in Arabidopsis thaliana, combinatorial interactions may result in many co-receptors with distinct auxin-sensing properties. We also demonstrate that the AFB5–Aux/IAA co-receptor selectively binds the auxinic herbicide picloram. This co-receptor system broadens the effective concentration range of the hormone and may contribute to the complexity of auxin response

Inter-domain Communication Mechanisms in an ABC Importer: A Molecular Dynamics Study of the MalFGK2E Complex

ATP-Binding Cassette transporters are ubiquitous membrane proteins that convert the energy from ATP-binding and hydrolysis into conformational changes of the transmembrane region to allow the translocation of substrates against their concentration gradient. Despite the large amount of structural and biochemical data available for this family, it is still not clear how the energy obtained from ATP hydrolysis in the ATPase domains is “transmitted” to the transmembrane domains. In this work, we focus our attention on the consequences of hydrolysis and inorganic phosphate exit in the maltose uptake system (MalFGK2E) from Escherichia coli. The prime goal is to identify and map the structural changes occurring during an ATP-hydrolytic cycle. For that, we use extensive molecular dynamics simulations to study three potential intermediate states (with 10 replicates each): an ATP-bound, an ADP plus inorganic phosphate-bound and an ADP-bound state. Our results show that the residues presenting major rearrangements are located in the A-loop, in the helical sub-domain, and in the “EAA motif” (especially in the “coupling helices” region). Additionally, in one of the simulations with ADP we were able to observe the opening of the NBD dimer accompanied by the dissociation of ADP from the ABC signature motif, but not from its corresponding P-loop motif. This work, together with several other MD studies, suggests a common communication mechanism both for importers and exporters, in which ATP-hydrolysis induces conformational changes in the helical sub-domain region, in turn transferred to the transmembrane domains via the “coupling helices”

Asymmetric Switching in a Homodimeric ABC Transporter: A Simulation Study

ABC transporters are a large family of membrane proteins involved in a variety of cellular processes, including multidrug and tumor resistance and ion channel regulation. Advances in the structural and functional understanding of ABC transporters have revealed that hydrolysis at the two canonical nucleotide-binding sites (NBSs) is co-operative and non-simultaneous. A conserved core architecture of bacterial and eukaryotic ABC exporters has been established, as exemplified by the crystal structure of the homodimeric multidrug exporter Sav1866. Currently, it is unclear how sequential ATP hydrolysis arises in a symmetric homodimeric transporter, since it implies at least transient asymmetry at the NBSs. We show by molecular dynamics simulation that the initially symmetric structure of Sav1866 readily undergoes asymmetric transitions at its NBSs in a pre-hydrolytic nucleotide configuration. MgATP-binding residues and a network of charged residues at the dimer interface are shown to form a sequence of putative molecular switches that allow ATP hydrolysis only at one NBS. We extend our findings to eukaryotic ABC exporters which often consist of two non-identical half-transporters, frequently with degeneracy substitutions at one of their two NBSs. Interestingly, many residues involved in asymmetric conformational switching in Sav1866 are substituted in degenerate eukaryotic NBS. This finding strengthens recent suggestions that the interplay of a consensus and a degenerate NBS in eukaroytic ABC proteins pre-determines the sequence of hydrolysis at the two NBSs

Activation of Akt by the Bacterial Inositol Phosphatase, SopB, is Wortmannin Insensitive

Salmonella enterica uses effector proteins translocated by a Type III Secretion System to invade epithelial cells. One of the invasion-associated effectors, SopB, is an inositol phosphatase that mediates sustained activation of the pro-survival kinase Akt in infected cells. Canonical activation of Akt involves membrane translocation and phosphorylation and is dependent on phosphatidyl inositide 3 kinase (PI3K). Here we have investigated these two distinct processes in Salmonella infected HeLa cells. Firstly, we found that SopB-dependent membrane translocation and phosphorylation of Akt are insensitive to the PI3K inhibitor wortmannin. Similarly, depletion of the PI3K regulatory subunits p85α and p85ß by RNAi had no inhibitory effect on SopB-dependent Akt phosphorylation. Nevertheless, SopB-dependent phosphorylation does depend on the Akt kinases, PDK1 and rictor-mTOR. Membrane translocation assays revealed a dependence on SopB for Akt recruitment to Salmonella ruffles and suggest that this is mediated by phosphoinositide (3,4) P2 rather than phosphoinositide (3,4,5) P3. Altogether these data demonstrate that Salmonella activates Akt via a wortmannin insensitive mechanism that is likely a class I PI3K-independent process that incorporates some essential elements of the canonical pathway

Molecular dynamics simulations and drug discovery

This review discusses the many roles atomistic computer simulations of macromolecular (for example, protein) receptors and their associated small-molecule ligands can play in drug discovery, including the identification of cryptic or allosteric binding sites, the enhancement of traditional virtual-screening methodologies, and the direct prediction of small-molecule binding energies. The limitations of current simulation methodologies, including the high computational costs and approximations of molecular forces required, are also discussed. With constant improvements in both computer power and algorithm design, the future of computer-aided drug design is promising; molecular dynamics simulations are likely to play an increasingly important role

Molecular dynamics simulations and membrane protein structure quality.

Despite a growing repertoire of membrane protein structures (currently approximately 120 unique structures), considerations of low resolution and crystallization in the absence of a lipid bilayer require the development of techniques to assess the global quality of membrane protein folds. This is also the case for assessment of, e.g. homology models of human membrane proteins based on structures of (distant) bacterial homologues. Molecular dynamics (MD) simulations may be used to help evaluate the quality of a membrane protein structure or model. We have used a structure of the bacterial ABC transporter MsbA which has the correct transmembrane helices but an incorrect handedness and topology of their packing to test simulation methods of quality assessment. An MD simulation of the MsbA model in a lipid bilayer is compared to a simulation of another bacterial ABC transporter, BtuCD. The latter structure has demonstrated good conformational stability in the same bilayer environment and over the same timescale (20 ns) as for the MsbA model simulation. A number of comparative analyses of the two simulations were performed to assess changes in the structural integrity of each protein. The results show a significant difference between the two simulations, chiefly due to the dramatic structural deformations of MsbA. We therefore propose that MD could become a useful quality control tool for membrane protein structural biology. In particular, it provides a way in which to explore the global conformational stability of a model membrane protein fold